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Abstract 

We consider the natural supersymmetry scenario in the framework of the R-parity conserving 
minimal supersymmetric standard model (called natural MSSM) and examine the observability 
of stop pair production at the LHC. We first scan the parameters of this scenario under various 
experimental constraints, including the SM-like Higgs boson mass, the indirect limits from precision 
electroweak data and B-decays. Then in the allowed parameter space we study the stop pair 
production at the LHC followed by the stop decay into a top quark plus a lightest neutralino or 
into a bottom quark plus a chargino. From detailed Monte Carlo simulations of the signals and 
backgrounds, we find the two decay modes are complementary to each other in probing the stop 
pair production, and the LHC with \/s = 14 TeV and 100 fbt luminosity is capable of discovering 
the stop predicted in natural MSSM up to 450 GeV. If no excess events were observed at the LHC, 


the 95% C.L. exclusion limits of the stop masses can reach around 537 GeV. 


PACS numbers: 14.80.Da,14.80.Ly,12.60.Jv 


I. INTRODUCTION 


Although the standard model (SM) has been successful in describing the existing exper- 
imental data, it is suffering from the hierarchy problem and new physics based on certain 
symmetry is widely expected to appear at TeV scale to stabilize the electroweak hierarchy 
against radiative corrections. This belief was further strengthened by the recent discovery 
of the Higgs boson at the Large Hadron Collider (LHC) with its mass determined around 
125 ee This mass value agrees well with the prediction of low energy supersymmetry 
(SUSY), which is so far the most promising new physics candidate. 

In SUSY, all known bosons and fermions have their supersymmetric partners, and the 
scalar top quarks (called stop t; with i = 1,2), as the top quark partners, can modify the 
property of the SM Higgs boson by exactly canceling out the dangerous quadratic divergence 
of the top quark loop. Obviously, the experimental determination of the stop properties is 
crucial to unravel the nature of supersymmetry in protecting the Higgs mass at the weak 
scale and thus solving the hierarchy problem. In fact, such activities have been carried out 
extensively at the hadronic colliders such as the LHC and the Tevatron "m but in contrast 
with the strong mass bounds (about 1 TeV) on the gluino and the first generation squarks 
, a relatively light stop (say about 300 GeV) can not be excluded. Nevertheless, it should 
be mentioned that the recently measured Higgs boson mass around 125 GeV may give some 
indications for the stop sector] B . In the popular MSSM with moderate tan @ and large 
ma, the Higgs mass is given by |6 
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m = MZ cos? 28 + TOG 73 aria ( - sin): (1) 
where v = 174 GeV, X, = A, — cot 8 and Msg is the average stop mass scale defined by 
Ms = mimi, This expression indicates that, for the heavier stop tj around 1 TeV as 
discussed above, the lighter stop tı must be heavier than about 200 GeV and |X;| must 
be larger than 1.5 TeV in order to push the Higgs mass up to 125 GeV B. About these 
constraints, one should keep in mind that they are independent of the decay modes of t, 
but on the other hand, they may be greatly weakened if there exists additional contribution 
to the Higgs mass A " 

On the theoretical side, there are good reasons to consider at least one stop significantly 
lighter than other squarks with a mass around several hundred GeV. Firstly, in some pop- 


ular grand unification models, supersymmetry breaking is usually assumed to transmit to 


the visible sector at a certain high energy scale, and then Yukawa contributions to the 
renormalization group evolution tend to reduce stop masses more than other squark masses. 
Secondly, the chiral mixing for certain flavor squarks is proportional to the mass of the cor- 
responding quark, and is therefore more sizable for stops. Such a mixing will further reduce 
the mass of the lighter stop. Thirdly, in the MSSM the minimization conditions of its Higgs 
potential imply " 


Mà (mån + Ea) (mi, -X)ta8 , 


LA FLUE Ee eee 2 
2 tan? 8 — 1 a (2) 


where mj; and mh, represent the weak scale soft SUSY breaking masses of the Higgs 
fields, u is the higgsino mass parameter, tan f = v,,/vg. X, and X arise from the radiative 


corrections to the Higgs potential and the dominant contribution to the X, is given by 


3Y? Tm. 
Du~ TE x mi, (v Qi 1) . (3) 


These two equations indicate that, if the individual terms on the right hand side of Eq. 


are comparable in magnitude so that the observed value of Mz is obtained without resorting 
to large cancelations, the natural values of u and mz, should be around 100 GeV and several 
hundred GeV respectively. Numerically, the requirement of X, < M2/2 (or X, < v?) leads 
to m; , upper bounded by about 500 GeV (or 1.5 TeV) [11]. Moreover, we note that a 
light stop is also phenomenologically needed by the electroweak baryogenesis and may 
be welcomed by the dark matter physics " In the MSSM, although the gluino contribute 
to mi at the one-loop level and to mj, at two-loop level, the corrections are proportional 
to m; and can be greatly enhanced by the large gluino mass|14]. In order to the keep the 
naturalness, we expect mg to be lighter than about 3 TeV for m; « 1.5 TeV. However, since 
the current results of searching for the supersymmetry indicate that a gluino with mass 
about 1 TeV can safely avoid the LHC constraints, we require 1 TeV < mg « 3 TeV in our 
calculation. 

Motivated by the theoretical preference and the results of the LHC search for SUSY, 
recently the natural MSSM scenario attracted broad attention "m which focuses on the 


following parameter space of the MSSM 14 ‘ 


e |u| S100 — 200 GeV and m;,, X 1 — 1.5 TeV as preferred by Eq. (2) and Eq.(3); 


e 1 TeV S m; S 3-4 TeV to escape the LHC constraint and at same time to avoid 
spoiling color symmetry; while the electroweak-ino masses may still be at sub- TeV 


scale; 


H 
e ma ~ mi, |? S |u| tan B as suggested by the relation m3 ~ 24? 4-m?,, -- m2, - X, 4- Y 


and Eq.(2); 


e m ,, ™M,, ~ 10 — 50 TeV to provide a decoupling solution to the SUSY flavor and 
CP problems. 


Since the stop are relatively light and sensitive to probing this scenario, there have been 
recently many theoretical studies on the collider signatures of the light directly produced 
stop in the R-parity conserving and violating MSSM|19]. For example, by using the top 
tagging technique, the sensitivity of stop searches were studied in the hadronic, semi-leptonic 
and di-leptonic channels|20 . In order to suppress the di-leptonic top backgrounds, the 
authors in Ref. explore some new kinematic observables developed from Mr, to improve 
the sensitivity of the stop searches. For the small mass splitting between stop and top, it is 
pointed that the rapidity difference and spin correlation of the daughter products from stops 
decay can be helpful to discover the signal|25]. When the stop mass is close to the lightest 
supersymmetric particle mass, the monojet signature from 741; + j production is expected 
to be useful in detecting the stopbd. If the stop mass is degenerate with the sum of the 
masses of its decay products, the searches based on missing transverse energy (E55 or Ér) 
have significant reach for stop masses above 175 GeV|27]. When the R-parity is violated, 
the decay modes of the stop will be very different from the ones in R-parity conserving 
MSSM, such as stop decaying to dilepton and trilepton final states 2]. We also noted that 
the constrains on the light stop in the natural SUSY have been discussed by using the results 
from sparticles searches at the LHC, and indicated that they were mild and can be safely 
avoided — P hd hd " 

In this work, we investigate the potential of the LHC in probing the lighter stop tı 
predicted by the natural MSSM with R-parity, which is based on some considerations: (i) 
Most of the studies of the stop searches have been carried out under some assumptions 
at the LHC in a model independent way or in simplified models. It will be meaningful 
to explore what might happen in a realistic model like MSSM under the current available 


experimental constraints; (ii) Due to the R-parity conservation, there will be sizeable missing 
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energy appearing in the sparticles productions and decays, which can be easily identified 
in the LHC data; (iii) One interesting phenomenological feature of the natural MSSM with 
R-parity is that both the lightest neutralino and the lighter chargino are Higgsino-like, and 
consequently tı always decays dominantly into tX? and by? with yt — XQ0W'*, which can 
greatly simplify the analysis of the f, detection at the LHC. 

For this purpose, we first scan the parameter space of the natural MSSM by considering 
various constraints in Sec. II. Then in Sec. III we discuss the observability of £4 through the 
direct stop pair production in the allowed parameter space by performing the Monte Carlo 
simulations for the channel pp — tit? — tX9tX? and the channel pp — fit; — btb. We 
will present their corresponding sensitivities for 8 TeV LHC and for 14 TeV LHC respectively. 


Finally in Sec. IV, we summarize the conclusions obtained in this work. 
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FIG. 1: Scatter plot of surviving samples in natural MSSM, projected on the plane of m;, versus 


|A;|/Mg. All the samples here predict the SM-like Higgs boson in the mass range 125 + 2 GeV. 


II. SCAN OVER THE PARAMETER SPACE 


Motivated by the natural MSSM, we scan the parameter space of the MSSM as follows: 


1 < tanp «60, 100 GeV < u < 200 GeV, 100 GeV < M» <1 TeV, 
|A| € 3 TeV, 100 GeV € (Mg,,Myg,) < 2 TeV, 90GeV<Ma<1TeV. (4) 
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FIG. 2: Same as Fig.1, but showing the relic density of the neutralino dark matter. 


For other unimportant parameters, we fix all the soft breaking parameters in the slepton 
sector and the first two generation sector at 10 TeV, and we assume A; = Ay, Mg, = Mp 
and M; : Mz = 1: 2 (inspired by the grand unification relation). In our scan, we consider 


following constraints: 


e We require the SM-like Higgs mass within the range 125 + 2 GeV. We use the code 
FeynHiggs2.8.6 ld to calculate the mass and the code HiggsBounds-3.8.0 to con- 


sider the experimental constraints on the Higgs sector of the natural MSSM. 


e Since the natural MSSM has important implications in B-physics 33. we use the 
code susy.flavor v2.0 to consider the constraints from the processes B — X. and 
Bsa) — i s 


e We consider indirect constraints from the precision electroweak observables such as pi, 


sin? 9! sp mw and R,. We use our own code for such calculation [4]. 


e We require the thermal relic density of the lightest neutralino (as the dark matter 
candidate) is below the WMAP value ld. We use the code MicrOmega v2.4 to 


calculate the density. 


After analyzing the surviving samples, we find they have two characters. One is that the 


Higgs mass of 125 + 2 GeV requires mz, > 220 GeV and there is a rather strong correlation 
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FIG. 3: Same as Fig[I] but showing the cross section of ttf production versus the stop mass. 


between mz and the ratio |A;|/Ms, as shown in FiglI] Here we further clarify that, if Mg, 
and Mg, are at sub- TeV scale, the minimum of mz, will be enhanced to about 300 GeV 1 
The other feature is that for most cases, the values of u are significantly smaller than Mj 
so that the lightest neutralino is higgsino-like. Fig.2 indicates that the surviving samples 
lie within two isolated regions. We checked that the lightest neutralino is bino-like in the 
left region and higgsino-like in the right region. Here the bino(higgsino)-like means it is 
still a mixed state but the dominant component is bino(higgsino). For the light neutralino 
dark matter(bino-like), the main annihilation channel is through exchanging Z boson. The 
annihilation cross section is roughly proportional to 1/ (4m? — mz). When the neutralino 
mass is about 50GeV—60GeV, the annihilation cross section may be very large, so the 
relic density will be less than 0.1. When the neutralino becomes heavier(60GeV —90GeV, 
neutralino is still bino-like), the annihilation cross section will drop. The relic density 
becomes large and even exceeds the WMAP value, and these samples are excluded. This 
is the reason for the gap between 60GeV —90GeV. When the neutralino goes on becoming 
heavy(>90 GeV), the dominant component of the neutralino will be higgsino. The coupling 
between neutralino and Higgs gets important and annihilation rate goes up, then the relic 
density drops. 

About the natural MSSM, we have two comments. One is that in this scenario the di- 
photon signal of the SM-like Higgs boson can hardly be enhanced to satisfy the requirement 
of the LHC data. This is because in the framework of the MSSM, there are only two cases 


which can enhance the di-photon rate, i.e. the small ae scenario a and the light 7 
scenario ld. " and in each case a large u is needed. In the Ref.|40], the authors pointed 
that the light stop with large couplings to Higgs boson in the SM+stop model can improve 
the SM fitting to the LHC and Tevatron data by enhancing l'(h — yy) and suppressing 
T(h — gg). However, we should note that it does not mean that the di-photon production 
rate can reach the measurement of the LHC in a concrete MSSM model since the reduction 
of l'(h — gg) is usually much stronger than the enhancement of I'(h — ^) for large values 
of X,(X, = A; — u tan 8)|41]. The other is that recently the ATLAS collaboration searched 
for the gluino-mediated stop pair production followed by the decay tı > bXi — blv, 
which set a lower bound mz, > 450 GeV |3]. This conclusion is not applicable to the our 
calculations since we take the gluino mass to be larger than 1TeV in the allowed parameters 


space of the natural SUSY. 


III. OBSERVABILITY OF STOP PAIR PRODUCTION AT THE LHC 


In this section we discuss the LHC potential of discovering the stop through the direct 
stop pair production in the natural MSSM at s = 8, 14 TeV. In FigBlwe show the pp > tt; 
production rate at the next leading order for the surviving samples. In getting this figure we 
used the package Prospino2.1 and the parton distribution function CTEQ6.6m with 
the renormalization scale ur and factorization scale wp setting to m;,. This figure indicates 
that the maximal values of the cross section reach 5.5 pb and 25.7 pb for the LHC with 
vs = 8 TeV and ys = 14 TeV respectively, and with the increase of the stop mass, the 
production rates drop rapidly. 

In Fig[4] we present various decay branching ratios of f; which are obtained by using 
the package SDECAY |44]. This figure indicates that for m;, € 320 GeV where the decay 
channel f, — tX? does not open up, tı decays into by} with a ratio of 100%, and as the stop 
becomes heavier, the branching ratios for tı — by? and tı — tX? may still be around 50%. 
In contrast, the branching ratios for t decays into t$a and by} are usually less than 20%. 

In the following we perform detailed Monte Carlo simulations to investigate the observ- 
ability of the direct stop pair production at the LHC. We concentrate on the semi-leptonic 
analysis with the b-tagging efficiency 40%, where the signal is consisted of four jets(at least 
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FIG. 4: Same as Fig] but showing the decay branching ratios of the stop. 


one b-jet), one lepton (e or u), and missing transverse energy. We first consider the process 
pp — bti — (Ext) (EXT) — (BET XT) jiki) or (bjii) (0€ vx). (5) 


From the ATLAS search for the signal tt + EZ"55 bi. we can see that the dominant SM back- 
ground after the E7"5* and Mr cuts is tt di-leptonic channel with one lost lepton and two 
additional jets from initial state radiation to fake the hadronic W. Another backgrounds in- 
clude tt semi-leptonic channel, tt di-leptonic channel with one 7 from top decay misidentified 
as a jet, W+jets and ttZ. Here we emphasize that the ttZ background becomes important 
for a heavy stop and should be considered in estimating the significance. In our calculation, 
we normalize the signal and the tt background to their NLO values , and simulate the 
signal and backgrounds by icone id interfaced with PYTHIA and Delphes to 
carry out the parton shower and fast detector simulation. We use the anti-k, algorithm 
with the distance parameter R = 0.4 to cluster jets and the MLM scheme kd to match our 
matrix element with parton shower. We checked that the shapes of the matched W 4-1,2,3 
partons are very similar, and for simplicity, we take W+2 jets samples in our calculations. 
In our calculations, since we employ the the variable MW defined in net {24}. we checked our 
results with theirs for the same parameters at ,/s = 7 TeV and found they were consistent 
with each other. 

In Fig] we show the distributions of Éy, the transverse mass Mr defined in and 
Mj for the backgrounds and our benchmark point m;, = 429 GeV and mgo = 110 GeV 
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FIG. 5: The distributions of @+4-jets+ with respect to Ér, Mr and ME for the signal pp > 
üt — txt and backgrounds. In the ME, distribution, we impose the cuts Py > 150 GeV and 
Mr > 150 GeV on the events and only display the di-leptonic tt, ttZ backgrounds and signal, 


where the events with wrong or no solution for M% are included in the last bin. 


TABLE I: The significance of stop pair production pp > t,t? > tx9tx9 for 100 fbt luminosity 


after imposing various cuts. Here we take m;, = 429 GeV and X = 110 GeV for illustration. 


Ey-cut (GeV)| Mr-cut (GeV)| MW -cut (GeV) | S/VB (8TeV) | S/V/B (14TeV) 


150 150 173 3.11 8.60 


with ys = 8 TeV (similar results are found for \/s = 14 TeV). This figure indicates that 
most events of W+jj and semi-leptonic tt backgrounds are characterized by Ær < 100 GeV 
and Mr < 100 GeV, and most events of the di-leptonic tt backgrounds are characterized 
by ME, < 170GeV, while a significant fraction of the signal may have larger Ér, Mr and 
MS. Fig[lalso indicates that the distributions of the t£Z background are quite similar to 
the signal and are difficult to be suppressed. Fortunately, the production rate of ttZ is much 
smaller than the one of t,t1. 

In Table I, we present the significance S/ VB of our benchmark point for 100 fb! lumi- 
nosity with \/s = 8 TeV and 14 TeV respectively by sequentially imposing the cuts on Ér, 


Mr and Ms. It can be seen that, for the given reference point, the cut Mr > 150 GeV can 
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greatly enhance the significance and MW, > 173 GeV further improves the significance by 
about 14% for J/s = 8 TeV and 9% for ,/s = 14 TeV to reach 3.11 and 8.60 respectively. 
Therefore, for our simulations in the allowed parameters space, we take the following 


events selection criteria: 


e One isolated electron or muon that passes the following requirements; 


— Electrons Ep > 25GeV and |n| < 2.47 without 1.37 < |n| < 1.52; 
— Muon: Er > 20GeV and |n| < 2.5; 


— Events are rejected if they contain a second lepton candidate with Pr > 15GeV; 
e Four or more reconstructed jets with Pr > 25GeV and |n| « 2.5. 


e Ér >150GeV, Mr >150GeV, MW >173GeV. 


where the basic cuts about pr and 7 on leptons and jets are from the ATLAS report}. In 
order to improve the signal sensitivity, we increase the values of ATLAS cut Ér from 100 
GeV to 150 GeV to further suppress the semi-leptonic tt background and use the new cut 


MW, >173 GeV to reduce the di-leptonic tt background in our calculations. 
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FIG. 6: The significance of stop pair production pp — tit} — tXQtXQ for the surviving samples in 


the natural MSSM. 


In Fig.6 we show the significance of the surviving samples with m; < 800 GeV. This 


figure indicates that the largest significance can be reached at m;, c 430 GeV where the 
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FIG. 7: Same as Fig[5] except that in plotting ME distribution we impose the cuts Æp > 100 


GeV and Mr > 100 GeV on the events. 
significance is about 1.5 for \/s = 8 TeV with 20 fb! luminosity and 8.5 for \/s = 14 TeV 
with 100 fb-! luminosity, and with the increase of the stop mass, the significance drop by 


are not as optimistic as those in bil because we have taken into account the branch 


one half for mz, ~ 500 GeV mainly due to the reduction of the production rate. Our results 
" 

ratio of 14 > tX9. Figl6lalso indicates that there are two branches for the significance in the 
mass region 320 GeV< m < 600 GeV. We checked that the upper branch corresponds to 
high branching ratio of tı — tX, which varies from 42.1% to 57.7% and results in a large 
signal rate, while the lower branch corresponds to a small ratio due to the competition of 
the decay mode t, — tX9 with tı — by} and t > tX}. 

Above analysis implies that, in order to fully explore the parameter space of the natural 


MSSM in stop detection, the decay mode tı — by{ should also be considered. So we next 


consider the process 

pp > titi > (bXx1)(bx1) > (0C vx1)(633X3) or (633X1) (0€ vX1). (6) 
Same as in Fig[5] we show the distributions of the three variables in Fig[7]for the benchmark 
point m; = 273.6 GeV, mx, = 163.5 GeV and myo = 156.3 GeV. Compared with the 
distribution in Fig[5] one can see that more signal events have lower values of # and 
low Mr, due to the relatively light t1. Fig[7] also indicates the MW, variable is helpless in 


suppressing the di-leptonic tt events and any cut on MJ) may hurt the signal greatly. 


In Figl8] we show the significance of the surviving samples for the process in Eq.(6). In 
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order to keep more signal events, here we relax the cuts of Ær and Mr used for the process 


in Eq.(5) as follows: 
Ep > 100 GeV, Mr > 100 GeV. (7) 


From this figure, one can learn that, due to the large stop pair production rate, the signifi- 
cance for mz, = 250 GeV may reach 7 for ys = 8 TeV with 20 fb"! luminosity and 64 for 
V/s = 14 TeV with 100 fbt luminosity, but for mz, = 400 GeV, the maximum value drops 
to 1.5 and 10 respectively. 
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FIG. 8: The significance of stop pair production pp — ftf — (bX1 (bx) for the surviving samples 
in the natural MSSM. 


Finally, we summary the significance of the direct stop pair production with the above 
two decay modes of tı for /s = 14GeV and 100 fb~! luminosity. The results are displayed 
in Fig[9] where only the maximal significance under each cut is shown. This figure indicates 
that, for mz, <400GeV, detecting the stop pair production through the chargino decay is 
more effective, while for 400GeV < m;, < 450GeV the neutralino decay is more effective. 
This figure also indicates that the LHC can discover tı predicated in nature MSSM up to 
450GeV. If no excess events were observed at the LHC, the 9596 C.L. exclusion limits of the 
stop masses can go up to around 537 GeV no matter what decay modes of the stop in the 


natural MSSM. 
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FIG. 9: The maximal significance in Eq.(5) or Eq.(6) as a function of the stop mass. The bullets 


(red) are obtained in Eq.(6) while the triangles (green) are obtained in Eq. (5). 


IV. CONCLUSION 


In this work we studied the direct stop pair production at the LHC in the natural MSSM. 
We scanned over the corresponding parameter space by considering various experimental 
constraints and then in the allowed parameter space we examined the observability of the 
direct stop pair production at the LHC through the semi-leptonic analysis. We focused on 


the following two channels 


pp — titi > (X1) (EX1) > (0C vx) (033X1) or (bj jX) (0€ vX1), 

pp > tt + (bxt) (ba) > (WL vX1)(633X1) or (b33X1) (0€ vX1). 
and performed detailed Monte Carlo simulations about the signals and backgrounds. We 
found that for m; < 400GeV the second channel is better while for 400GeV < m; < 
450GeV the first channel is better. We also found that the LHC with ys = 14TeV and 
100 fb-! luminosity is capable of discovering tı predicated in nature MSSM up to 450GeV. 
If no excess events were observed at the LHC, the 9596 C.L. exclusion limits of the stop 


masses can reach around 537 GeV in the natural MSSM. 
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FIG. 10: The ATLAS observed and our expected exclusion limit in the plane of mgo versus mg, 


assuming BR(t, > tx?) =100%. 
Note Added 


Very recently, the ATLAS collaboration reported the result of the direct searching for the 
stop pair production base on 4.7fb~! of data|b1]. We validated our simulation by reproducing 


the ATLAS exclusion limit according to the assumptions and cuts in the report as follows: 


e One isolated electron or muon passing ‘tight’ selection criteria; 


— Electrons Ep > 25GeV and |n| < 2.47; 


— Muon: Er > 20GeV and |n| < 2.4. 


e Four or more jets with |n| < 2.5 and Pr >80, 60, 40 and 25 GeV, and at least one jet 
to be identified as a b-jet; 


e Admin > 0.8,where Admin is the minimum azimuthal separation between the two 


highest Pr jets and the missing transverse momentum direction; 
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The jet-jet pair having invariant mass > 60 GeV and the smallest AR is selected to 
form the hadronically decaying W boson. The mass mj; is reconstructed including 
a third jet closest in AR to the hadronic W boson momentum vector and 130 GeV 


< mj; < 205 GeV is required; 


ERs 2150 GeV, ET55// Hy >7 GeV"? and mr > 120 GeV; 


Events are rejected if they contain additional leptons passing looser selection criteria. 


Here we treat the looser selection criteria as Pr > 15 GeV. 


The branch ratio of tı — tX? is assumed to be 100%. 


In Fig.10, we display the expected exclusion limit from our simulation. Considering the 
differences between the fast simulation and full detector simulation, we can see that our 
result is consistent with the ATLAS exclusion limit within the reasonable error range. We 
also expect our result can be improved by the simultaneous fits method used by ATLAS 
for five signal regions and three control regions, however, which is beyond the scope of our 
simulation. It should be noted that the above stop masses limits can be avoided in our 


study, since the stop decays with a mixture of the branching ratios. 
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